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Abstract 

We investigate periodic straight-line orbits (SLO) in Hamiltonian force fields using both 
direct and inverse methods. A general theorem is proven for natural Hamiltonians quadratic 
in the momenta in arbitrary dimension and specialized to two and three dimension. Next we 
specialize to homogeneous potentials and their superpositions, including the familiar Hcnon- 
Heiles problem. It is shown that SLO's can exist for arbitrary finite superpositions of iV- 
forms. The results are applied to a family of generalized Henon-Heiles potentials having discrete 
rotational symmetry. SLO's are also found for superpositions of these potentials. 

1 Introduction 

The connection between the geometry of trajectories and the force fields that generate them has 
long been of interest to physicists and mathematicians alike |Sze67| IBoz95| IVDM91] . In the copious 
literature on the subject one may distinguish between direct and inverse approaches. The direct 
problem asks: what are closed form solutions for particular orbits in a given potential? By contrast, 
the inverse problem poses the question: which force field will produce a set of orbits with a given 
shape? For example, suppose that the "natural" Hamiltonian, H : M^" M, 

H{q,p) = ^\p\^ + V{q) , (1) 

has an orbit with energy E that lies on a given surface in the configuration space, 

F{q,c) = f{q)-c = . 



where / : M. As was first shown by Szebehely [Sze74| for the planar case, and generalized by 

Puel |Pue84] to n-dimensions, the potential must satisfy the equation 

^^W.2(E-^,V^(^) (2) 

on F = 0. If the orbit is restricted to a curve, given by the intersection of n — 1 surfaces Fi = 
fiio) ~ Ci, then the potential must satisfy Q for each i = 1, . . .n — 1. The potentials obtained 
in this way have an n — 1 parameter family of orbits lying on the intersection of the sets Fi = 0; 
these equations have been studied in a number of papers and many three dimensional examples 
have been obtained jPueSSl IBK04llBK05] . 

In contrast to Szebehely's problem, we consider the problem of finding a potential with a single 
orbit of a given shape. In lieu of specifying this curve as the intersection of surfaces, we find it 
more convenient to represent it parametrically through Q : M ^ M", as 

q{t) = Q{x{t)) , (3) 

where x{t) G M now represents the temporal dynamics. Note that x — Q{x) = may be regarded 
as a configurational invariant |Hal83| valid for one particular value of the energy E. 

Substituting this into the equations of motion q = —W gives Q'x + Q"x'^ = —W. Since 
bp = IQ'px^ = 2{E — V), we can eliminate the first derivative to obtain 



/O'" V7T/ ~ ^ n" 

Qx=-VV Q . 



2{E-V) 

The implication is that the vector on the right must be in the Q' direction, in other words that the 
projection onto the plane orthogonal to the tangent vector Q' must be zero. The projection matrix 
orthogonal to the vector Q' is 

(note that P"^ = P and PQ' = 0). We then obtain 

^>(vv+'^^Q")=o (4) 

along q = Q{x). If V satisfies Q then the dynamics reduces to the scalar system 

- = -%i ■ v^(Q) - 2%:r^{E - vm . (5) 



Thus the inverse problem reduces to finding a potential that satisfies Q. In general this seems to 
be a hard problem that we will leave to a later paper. 

The simplest geometry for an orbit is a straight line in the configuration space: 

q{t) = Q{x{t)) = sx{t) + go , (6) 

for a constant "slope vector" s, i.e., a straight line orbit (SLO). In this case, Q simplifies con- 
siderably since Q" = 0, and the requirement on the potential is simply that its gradient must be 
parallel to Q' = s, or specifically 

VV\sx+qo = A(x)s , (7) 
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for some scalar function A. In this case the dynamical equation (|5]) reduces to 

X = -\{x) . (8) 

Since ([7]) is independent of the energy, the straight line orbits that we find automatically come in 
one-parameter families, parametrized by E. 

Our investigations possess some of the aspects of both direct and inverse methods. Indeed, for 
a given potential V , we can solve ([7| to determine the allowed values of s, if any. We will call 
the set of admissible slopes, the slope spectrum of V . Alternatively, we can fix s and solve the 
eigenvalue-like problem ([T]) for the potential; the general solution of ([7| will be obtained in ^ 
Subsequently, we will specialize to the two and three degree-of-freedom cases, giving a number of 
examples. For two degrees of freedom, the general solution to ([7| involves two arbitrary functions, 
and for three, it involves four functions. Another class of examples, superpositions of homogeneous 
potentials, is treated in ^ An example of this case is the famous Henon-Heiles system which has 
three families of SLOs. By expressing this case in polar coordinates we also obtain SLOs for a 
family of hyper-Henon-Heiles potentials. 

In all cases we are motivated by physical applications and decline interest in discovering exotic 
potentials which will never be found in a physical problem. 

The ideas here should be contrasted with the notions of central configuration and choreography 
in celestial mechanics. A central configuration is a solution in which q = —X{q)q, for some scalar 
function A : M |! Moe90] . When the masses are equal, the potential must satisfy W = X{q)q, 

instead of Q or ([7]). A choreography is a solution of an A^-body system in which each body has 
identical configuration space, and each follows the same curve, but with a phase shift |CGMS02] . 
In the standard gravitational problem, the configuration space is M'^^ for bodies in and each 
body has a configuration orbit that lies on a curve qi{t) = C{x{t) + 4>i) for a curve C : M — > M^. 



2 Straight Line Orbits 

As in the introduction, we consider a Hamiltonian system on M^" with coordinates {q,p). A straight 
line orbit has the form ([6]) with intercept ^Oi slope vector s G M" and scalar dynamical function 
x(t). For any particular SLO we can, without loss of generality, choose coordinates so that go = 0; 
consequently, we will look only for orbits that go through the origin. Moreover, since the equation 
is homogeneous in s, we can choose the slope vector so that s| = 1. 

In particular consider the natural Hamiltonian system with potential V G C^(M",M). If H 
admits an SLO with slope s, then 

p = sx 
sx = -W{q)\g=sx ■ 

Thus an SLO exists only when W^sx) is parallel to s for all x G M. We call the set of admis- 
sible slopes the slope spectrum, Sl{H) of the Hamiltonian H. The slope spectrum for ([T| is thus 
determined by a nonlinear eigenvector-like equation 

Sl{H) = {s G M" : VV{sx) = X{x)s , A : M ^ M}. (9) 

The general form of a potential admitting an orbit with a given slope can be easily determined: 
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Theorem 1. The Hamiltonian system ([Ip has a family of straight line orbits q{t) = sx{t), \s\ = 1, 
only if the potential has the form 



where W - 
In this case x : 



V{q) = U{s^q) + \q^W{q)q , (10) 

is any symmetric matrix function that has a zero eigenvector s, W{q)s = 0. 
is any solution of the one- dimensional ODE 

-U'{x) . 



X 



Proof. It is not hard to see that (10) satisfies ([9|). To show this is the general form, we choose a 
new basis aligned with s. Let R = (s, M) be an n x n orthogonal matrix, so that the columns of 
the n X (n — 1) matrix M are orthonormal and orthogonal to s: M'^s = 0. The slope requirement 
^ then becomes the system of n — 1 equations: M'^W{sx) = 0. 

Defining new coordinates by g = i?^ = sx + My, where ^ = {x,y) and y G M""-*^, then the 
potential in the new coordinates is U{x,y) = V{sx + My). Noting that Vq = R^^V^ = RV^, the 
slope equation becomes 



M^RV^^^y)Uix,y) 



y=o 



. 



But since R is orthogonal. 



so M'^R = (0, /) (where this I has size n — 1) and the slope equation reduces simply to 



dyU{x,y) 



. 



:i2) 



This is just the requirement that U has a zero derivative with respect to the n — 1 variables y when 
they vanish. It has general solution 

U{x,y) = U{x) + ly^W{x,y)y 

where is a smooth, symmetric (n — 1) x (n — 1) matrix function. In terms of the original 
coordinates, note that y = M^q, and x = s'^q so that 



V{q) = U{R^ 



U{s'^q) + \q^MWM^q . 



Note that the matrix W = MWM is symmetric and that its rank is no more than n—1; it has 
a zero eigenvector s. Thus we have the general solution (10). The equation of motion on the SLO 
immediately reduces to the one-degree-of- freedom system (11). □ 



If V has a local minimum at the origin, then (11) will have some bounded solutions. However, 
this does not imply that the resulting orbits are stable when thought of as orbits of the full system. 
We turn next to some examples for two and three degrees of freedom. 
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2.1 Two-degree-of-freedom natural flows 



Here we give an explicit form for ( 10 ) for the case of two degrees of freedom assuming a straight 
Hne orbit of the form q(t) = {x(t), ax{t)). In this case, the requirement on the potential reduces to 
the single equation 

Sm = |a : - V{x, ax) = o} , (13) 

which could easily be solved directly. However, it is also a simple application of Th. [l} 

Corollary 2. The natural Hamiltonian |Ip with {q,p) = {x,y,Px,Py) has a straight line orbit of 
the form y{t) = ax{t) only if 

V{x, y) = F{x + ay) + {ax - y)^G{x, y) , (14) 

where G is continuous at y = ax. In this case x obeys the ODE x = —F'((l + a'^)x). 

Proof. For this case s oc (1, a) and the coordinate transformation is 

n / 1, ^\ 1 /I —a 

R = {s,M) 



Vl + a^ \a 1 



and the 2x2 matrix W m ( 10 1 has s as a zero eigenvector when it is proportional to MM"'", giving 

W{q)=2F{x,y) 



o? —a 
—a 1 



Thus (10) becomes 

V{q) = JJ{s^ q) {ax - yfF{x,y). 



Scaling the argument of U gives ( 14 ) □ 
For example 

V{x, y) = cos(x + y) + (x — 7/)^ sin x cos y 
has orbits y(t) = x(t) with x obeying the pendulum equation 

X = sin 2x. 

It is also easy to construct potentials which have multiple straight line orbits. For example: 

V{x,y) = {x- yfF{x,y) , 

has X = y as an orbit. We now may replace i*" by a function that has other straight lines; for 
example, 

y(x, y) = (x - yfx^ix + yf 
has y = X, X = and y = — x as orbits. 
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2.2 Three-dimensional natural flows 



Consider the three degree of freedom case of ([T]) and a straight hne orbit of the form L 
{{x{t)^ax{t)^ I3z{t))}. The slope spectrum is then given by 



Sl{H) 



a ay 



pdz 



(15) 



This could be solved directly, but it is also easy to directly apply Th. [T} 

Corollary 3. The natural Hamiltonian |ip with {q,p) = {x,y, z,px,Py,Pz) has a family of straight 
line orbits q{t) = {x{t), ax{t), Px{t)) only if functions A(x,y,z), B{x,y,z), C{x,y,z) and U{x) 
exist such that 

V{x, y, z) = U{x + ay + f3z) + (ax - yfA + (/3x - zfB + (/3y - az^C , 
where A is continuous at y = ax, B at z = fjx and C at (3y = az. 



W 



Proof. Here we set s oc (l,a,/3). Requiring that the matrix W in (10) has s as a zero eigenvector 
yields the form 

^a'^A + 13'^B -aA -/SB 
-aA A + P^C -apC 
-pB -ape B + a'^Cj 

This results in the quadratic form 

q^Wq = {ax - yf A + (/3x - zfB + {Py - azfC 



which immediately gives the result. We will give several examples in Sec. 4. 



□ 



3 Homogeneous Potentials 



Homogeneous potentials are often encountered practice, and we are therefore motivated to develop 
a method specifically tailored for this class. In particular suppose that H has the form 



H 



\\p\' + Y.U^'\q) , 

k=2 



(16) 



where each term U^^\aq) = a'^U^^\q) is homogeneous with degree k. We shall concentrate on the 
two and three degree of freedom cases. Examples include the Henon-Heiles system |HH64j and its 
generalizations |Hal83] . where the potentials are polynomial. 

Note that both (13) and (15) reduce to individual equations for each homogeneous term in U: 



Lemma 4. The slope spectrum for (16) is the intersection of the slope spectra for the homogeneous 
potential system H^^^ = + U^^^{q). For any s £ Sl{H), with \s\ = 1, the orbit q(t) = sx{t) 

satisfies 



N 



sjx 



(17) 



k=2 
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Proof. The requirement of (|9|) becomes 



N 



^(VC/W(s)-A,„i5)x*=-i = 0, 



k=2 



where A(x) = A^x^ is also a polynomial in x. Each term in the sum is a homogeneous 

polynomial in the scalar function x{t) of degree k — 1, and, unless x{t) is constant, these terms 
must vanish individually since different powers of a non-constant function are linearly independent. 
This gives the individual "eigenvalue" problems 



(18) 



Thus 



N 



Sl{H) = Pi 5/(F(^)) . 



k=2 



As in the general case, for each s £ Sl{H), x must solve (g. By Ak-i = s ■ VU^''\s), 
which by homogeneity becomes Afc„i = kU^^\s). Thus te| reduces to (17). □ 



Specializing now to the case of polynomials, we consider first the 2D case where 

k 

U^^\x,y)=^anx^y^-^ . (19) 



n=0 



If we look for an orbit with slope s = (1, a), the slope spectrum requirement (13 1 reduces to 



N 



Q^{a) = ^[na„a - {N - n + l)an-i] 



a^-" = . 



(20) 



n=l 



This equation can be thought of in two ways. For a given a, it can be viewed as a single linear 
restriction on the coefficients. Alternatively, for a given set of coefficients, (20) becomes a single 
polynomial equation in a whose real zeros determine the slope spectrum. 
In particular the quadratic case 

Q^'^\a) = a\o? + 2(a2 — ao)a — a\ 

always has two real zeros since its discriminant 

A = (a2 - ao)^ + a\ 

is nonnegative. There is one special case, in which any a is in the slope spectrum: Q^^^ is identically 
zero when ai = and 02 = oq, which corresponds to the harmonic oscillator case 

^suo = \{x^^y') ■ 

The cubic case reduces to 

(5^^^(a) = a\o? + (2a2 — 3ao)a^ + (803 — 2a\)a — 02. 
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This, of course, has at least one real zero, so the slope spectrum is always nonempty. 

It is also easy to find examples of homogeneous, but non-polynomial potentials with SLO's. For 
example, for the degree-two potential 



we obtain 



9x^ + 8x^y + 



g(2)(«)=4fc«)(5«-l)(l + «) 



which has three real zeros giving the slope spectrum 

Sl{H) = {2,-1,1} 
Now consider the three degree of freedom case, with 



= Yl anmgx^y"'z'^. (21) 

m+n+q=N 



The requirement (15) reduces to the two equations 

Q[''\a,l3) = Y,anmg{ma"'-' - na"'+^)(3'i = 



(22) 



Again these equations can be viewed in two ways. For a given pair (a, /5) they are a set of si- 
multaneous linear equations for the coefficients amnq- Alternatively, for given coefficients, the two 
polynomials many have a set of simultaneous solutions that give the slope spectrum. These solutions 
may be found by taking the resultant of Qi and Q2- 



4 Examples 

4.1 The Henon-Heiles System 

The well-studied Henon-Heiles system is a natural Hamiltonian with potential |HH64| 

y(x,y) = i(x2 + y2)+^2y_iy3_ (23) 

The standard {y,Py) section with a; = and px > for total energy E = H = ^ is shown in Fig. [ij 
To find straight line orbits for the potential (23), we can apply Lem.|4j Since the quadratic part 
is symmetric, it has SLOs for any s. The slope spectrum ^ for the cubic part gives the eigenvector 
equation 

' = A 

which reduces to the equations si = and S2 = —A, or sf = 3^2 and 2s2 = A. Thus the slope 
spectrum for the Henon-Heiles Hamiltonian is 

Sl{H) = {{0,1), iV3,l), {-V3,l)} . 
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Figure 1: Poincare section for Henon-Heiles potential, with E = 1/9, 



These three straight hne orbits are shown projected onto configuration space in Fig. [2ja). Once 
we know the slopes, we can show that (23) is of the form (14) by first solving for F using 



F ((1 + a^)x) = V{x, ax), and then solving for G. For example for q = we find 

F{x) = ^x2(x + ^/3) , 

o 

3 

G{x,y) = g(i - - y) ■ 

The SLO orbits with a = ±3"^/^ appear as fixed points on the y— axis in the section of Fig. [1} 
These orbits cross the section x = at y = and since py = apx, this implies that the momentum 
on the section is 

For E = |, the SLOs are elliptic and form the centers of the two stable islands at py = ±\ in 
Fig. [T] These fixed points are stable up to ii^ ~ 0.14, at which point a pitchfork bifurcation occurs 
generating a pair of stable periodic orbits that are no longer SLOs; one of these orbits is shown in 
Fig.|2|b). The vertical SLO, with x = lies in the section and corresponds to its boundary, namely 
the contour of 

E = H{0,y) = l{pl + y^)-yK 

Similar results have been obtained previously by Antonov and Timoshkova |AT93] and van der 
Merwe |VDM91j . 
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Figure 2: (a) The three SLOs for the Henon-Heiles system for E ^ 1/8, (b) stable bifurcated orbit at E ^ 0.14. 



4.2 A Quartic Potential 

As is well known, the potential for the Henon-Heiles system can be written in polar coordinates as 



■ sin 3(j). This suggests the quartic analogue 



■ sin 4<j) 



(25) 



For this potential (20) becomes the biquadratic, 

Q(4)(a) =a^-Qa^ + 1 = 

with roots = 3 it 2-v/2 which implies that a = zttanvr/S and ibtan37r/8. The resulting four 
SLOs are shown in Fig. ^a) for E = ^. 

The section x = 0, Px > is shown in Fig.|4j for E = ^ the four SLOs correspond to elliptic fixed 
points on the y-axis at the momenta determined by (24), namely, pyo ~ ±0.191 and ±0.462. As the 



energy increases, bifurcations occur, resulting in a changing number of fixed points orbits on y-axis. 
These orbits are ephemeral, coming and going with changing E. For example at E = the section. 



shown in Fig. |4|b), there are at least six elliptic fixed points on the y-axis. In Fig. [3][^b)-(c) two 
of the additional period-one orbits, located at pyQ ~ —0.213 and —0.2518, are shown. Although 
these are clearly not SLOs, the first one is remarkably linear near the origin and self-retracing. 
Fig. [5] depicts a bifurcation sequence that generates additional period-one orbits. For E = 0.21 
there are two SLOs in the lower half-plane; for E = 0.22, one SLO has destabilized by a subcritical 
pitchfork bifurcation, spawning two stable non-SLOs. At E = 0.23 this SLO has restabilized via a 
supercritical pitchfork bifurcation. There are now a total of 4 period-one orbits in the lower half 
plane, of which two are SLOs. Of course, the number of SLO's is constant. 
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4.3 Polar Coordinates 

The polar coordinate construction in the previous example suggests the following 
Lemma 5. The two degree- of- freedom Hamiltonian 



1 2 



n2 



H=-pi + ^ + V{r,4>), (26) 
where V is smooth, has an SLO (r(t),(/)o) if 

V = R{r)+P{<l,)S(r,(t>) , (27) 
where P'{(j)Q) = and R{r) and 5(r, (p) are arbitrary smooth functions. 



Proof. From the equations of motion for (26) in polar coordinates (j) = when 

dV{r,cj)) ^ 

d4> ''^=^" 



with general solution (27). □ 
As an example, consider the /lyper-Henon-Heiles family of potentials 

[/W(r,(/)) = — sinA^0 (28) 

which has Cn symmetry for any positive integer A^. Including the harmonic potential gives the 
corresponding Hamiltonian 

1 1 
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Figure 4: Poincare section for quartic potential (251, with (a) E ^ 1/8, (b) S = 1/4 



where 



h. Straightline orbits correspond to zeros of 



V4> 



dH 



-r^ cos N(p 



so that Ncj) = n7r/2, with n a positive odd integer. For the Henon-Heiles system {N = 3), (pi = gvr, 
jvr, and gTr, in agreement with Fig. [2]|^a). For the quartic system {N = 4), 4>i = ^ir, ^tt, gTr, and 
gTr in agreement with Fig. [s^a). The motion on an SLO is given by 



f. + ^ + (_l){"-i)/V^-i = 0. 

Note that the orbits with n = 3 mod 4 will be unbounded if their energy exceeds the threshold 
E = ^^jf^, but the orbits with n = 1 mod 4 are bounded for all positive energy values. 



Finally, consider the superposition of two hyper-Henon-Heiles potentials 



jjiN) _^ ij{M) ^ !_ sin AT^ + !_ sin M(t>. 



(29) 



By Lem. |4j this potential has SLOs at the common slopes of the two homogeneous potentials, or 
when 

nvr mvr 
~ 2iV ~ 2M ' 

for some positive, odd integers m and n. Hence, SLO's occur whenever nM = niN, an interesting 
little problem in Diophantine analysis. Since (n, m) are odd it is clear that and M must both 



be even or both odd. Thus, a superposition of the Henon-Heiles potential (23) and the quartic (25 ) 
has no SLOs. If N and M are both odd, a simple family of solutions occurs when n = kM, m = kN 
for any natural number k. For example, for {N,M) = (3,5), an SLO occurs for {n,m) = (5,3). 
Note that if N and M have any common factors then these can be removed from the homogeneous 
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Figure 5: Bifurcation sequence near py — —0.2 for quartic potential (25 I. (a) E = 0.21 single stable SLO, (b) 
E — 0.22 SLO destabilized via subcritical pitchfork bifurcation, (c) E = 0.23 SLO restabilized via supercritical 
pitchfork bifurcation. 



equation, and once they are removed the remaining factors of -A'^ and M must both be odd (if they 
were both even, then another factor of 2 can be removed). Thus when .A'^ and M are both even 
then they must have the same power of two in their prime factorization. For example, if = 2^, 
then we must have M = 2^(2j + 1) for some integer j. Thus the first common SLOs occur when 
M = 12, for example, with (n,m) = {k,3k). 

4.4 Three Dimensional Examples 



The direct problem (22 ) is easily solved for any given potential. For example, the 3D Henon-Heiles - 
like model 

1/9 O 9^ / 9 9^ 1 



V{x, y, z) = ^(x2 + y2 + z^) + (x^ + z^)y - -{y^ + z^) , 



has two SLOs: 



^W = {(^,0),(-^,0)} 



Similarly, the potential 



V{x, y, z) = -(x^ + + z^) + + y^ + z^ + xyz 

has five SLOs: 

5/(i7) = {(0,0),(i,l),(l,i),(l,l),(2,2)}. 

Since for any degree, ( |22[ ) can be viewed as just two equations for the coefficients of F as a 
function of a and /?, there are many solutions of the inverse problem. For example, the cubic 
potential 

1 _ 2/32 _ 3 , ^ B 



3 I 2 I M 2 

-X +yx -\ X z . 

3a a 



has an SLO y = ax, and z = (3x. When a = (3 = 1 this reduces to 

^(3) ^ _^3 ^^y^ ^^^2 _ 



(30) 
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This potential has an additional SLO, that can be found from (22), so that its slope spectrum is 

5/(F) = {(1,1), (1,-1)}. 
A contour plot of U for E = 0.1, including the SLO's is shown in Fig. |6j 




Figure 6: Contour plot of the cubic potential (30 1 for E = 0.1, showing straight line orbits 



A quartic solution of (22) is 



So, we can get a nice example, we compact energy surfaces, and presumably chaotic orbits by 
putting these together 



H 



Since the quartic terms dominate for large coordinates, they bound the motion. 

As in the 2D problem one can readily incorporate multiple SLOs in the inverse problem. For 
example, from Cor. [3] the potential 

V{x, y, z) = x^Fi{x, y, z) + Gi(y, z) = /F2(x, y, z) + G2{x, z) = z^F2{x, y, z) + G3{x, y). 

has three SLOs along the coordinate axes. That is, we want the equations of motion for each 
variable to be of the form x = xK{x, y, z), so that x = is a solution. An example is 



V{x, y, z) = ax^ + I3y^ + 72;^ + ax"^ + by'^ + cz'^ + dx^y^ + ex^z^ + fy'^z^ 
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5 Discussion 



We have determined very general conditions for SLOs for natural potentials in arbitrary dimension. 
Using these results one can either construct potentials with a given SLO or test a given potential for 
SLOs. The general solution for two degrees of freedom involves two arbitrary functions, for three 
degrees of freedom, four. Superpositions of potentials having SLOs are also easily constructed. 
The special case of homogeneous functions occurs rather frequently and as examples we studied 
the Henon-Heiles system and a family of hyper-Henon-Heiles systems, in polar coordinates. Several 
two- and three-dimensional examples have been analyzed. 

It may be possible to apply similar methodology to more general problems, e.g., to find all 
potentials with quadratic orbits and to generalize the form of the Hamiltonian to include a mass 
matrix. It would be interesting to learn whether similar behavior also occurs in non-Hamiltonian, 
reversible systems. 
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